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Abstract 

This paper explores an optical interpretation of the gravitational redshift effect, and 
shows how the deflection of light and the radar echo delay can be calculated from the point 
of view of optics. A formal equivalence between general relativity and geometrical optics, 
apparently unrelated areas of physics, is discussed. The development of an optical analogy 
suggests introducing an optical metric, an optical version of the Schwarzschild metric. 



1 Introduction 

The general theory of relativity has been accepted as the most satisfactory theory of gravita- 
tion. This acceptance rests partly on its conceptual and structural elegance, and partly on its 
agreement with experimental observation. But the theory has been essentially mathematical 
in character, being concerned with the consequences of a "geometrization" of the space-time 
manifold. Thus the application of the theory involves the use of special mathematical methods 
which, although relevant to optics in many cases, may easily be considered separately from it. 
In attempting to deduce their optical nature, one may pass from the mathematical language to 
the physical language, and see how they can be reconciled with each other. 

Four classic tests are usually cited as experimental verifications of the general theory of 
relativity: the gravitational redshift of spectral lines, the deflection of light by the Sun, the 
precession of the perihelion of the orbit of the planet Mercury, and the time delay of radar 
echoes passing close to the Sun. Three of these tests examine the influence of the gravitational 
potential on the propagation of light. Only the planetary orbit precession investigates the 
motion of a particle of finite mass in the gravitational field of the Sun. Because these are 
optical phenomena, one may raise a question as to whether the three classic tests can also be 
inferred from the point of view of optics. 

It is actually possible to predict these tests in a valid manner on the basis of optics. This 
paper explores an optical version of gravitational redshift, and shows how the deflection of light 
and the radar echo delay can be calculated therefrom. The interpretation proposed in this 
paper ascribes the redshift effect to an optical phenomenon related to the speed of light in a 
non-uniform medium. It will be of particular interest to note that the deflection of light and the 
radar echo delay can also be correctly derived from the equation for rays in geometrical optics 
without using the geodesic equations or the field equations of general relativity. As a matter 
of fact, there is a formal equivalence of expressions between general relativity and geometrical 
optics, apparently unrelated areas of physics. The development of an optical analogy suggests 
introducing an optical metric, an optical version of the Schwarzschild metric. 
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2 Redshift of Spectral Lines 

In 1911, Einstein' 1 ) predicted a change in the frequency of spectral lines with gravitational 
potential, generally referred to as the gravitational redshift. The argument, as Einstein explains, 
was as follows: 

"Let the two points S and S' be placed at rest, a distance h apart along the lines 
of force in a uniform gravitational field of acceleration g. In accordance with the 
principle of equivalence, we are able, in place of the system K in this gravitational 
field, to set the gravitational-free system K' which is accelerated with — g. Consider 
the process of propagation of radiation from S to S' from a system Kq, which is 
to be free from acceleration. At the moment when the radiation of frequency / is 
emitted from S to S', let the velocity of K' relative to Kq be zero. The radiation 
will arrive at S' when the time h/c has elapsed to a first approximation. But at 
this moment the velocity of S' relative to Kq is v = gh/c. Therefore, by Doppler's 
principle, the radiation arriving at S' does not possess the frequency / but a greater 
frequency /' which is related to / to a first approximation by the equation 

By the equivalence of K and K' , we may replace gh by the gravitational potential 
if the same process takes place in the system K." 

Astronomical observations, though somewhat ambiguous, have tended to confirm this effect. 
Since it does not seem possible to predict this interesting effect without using the general 
theory of relativity, the gravitational redshift is now recognized as resulting from the principle 
of equivalence. Contrary to the current recognition, however, it would always seem possible 
to find a natural expression for the gravitational redshift from the optical point of view. If 
we seek to find an optical nature of the gravitational redshift, we must turn our attention to 
the velocity v = gh/c, not to the principle of equivalence. At least phenomenologically, the 
effect would appear to be due to this change of velocity which the radiation experiences during 
the propagation along the lines of force of the gravitational field. Looking for a Newtonian 
mechanical interpretation of it, one finds without difficulty that it is equal to the velocity 
difference due to the medium or fluid which the radiation experiences during the propagation 
between the places at different gravitational potential. 

Let p be the density and 4> be the gravitational potential with g = — V</>. A pressure 
uniform throughout a fluid mass produces no effect on the motion. The time rate of change of 
the momentum of fluid is equal to and opposite to the pressure gradient force in the medium. If 
we calculate the velocity difference due to fluid with differing gravitational potential according 
to 

dv „ , 

— = -g or V0, (2) 

we obtain the same velocity as that in (1). This leads to a simple physical interpretation of 
it: the redshift effect is attributed to the relative velocity change due to the fluid or medium 
by which light is affected during propagation. In such an interpretation, the redshift appears 
to be an optical phenomenon related to the velocity of propagation of light in a non-uniform 
medium. This leads us to consider the redshift effect as being purely optical in origin. There is 
no difference in form between the optical interpretation and Einstein's argument here quoted. 
In principle, there is thus no objection to interpreting the redshift result from the optical point 
of view. Although Einstein further explained it as due to the time dilatation between clocks in 
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different gravitational potentials, what he had found was phenomenologically no more than the 
fact that the velocity of light is altered, linearly to a first approximation, by the medium as a 
result of pressure gradient force. 

A difference of interpretation already existed at the first time of observation. Jewell in 1897 
and particularly Fabry and Boisson^ 2 ) in 1909 found displacements of solar spectral lines toward 
the red end of the spectrum, and ascribed them to an effect of pressure in the absorbing layer. 
However, Einstein's theory^ 3 ) of general relativity in 1916 established in most physicist's minds 
the interpretation of redshift as a manifestation of time dilatation in a gravitational potential, 
and this rather unusual explanation has survived until the present. Like most of the relativistic 
explanations, it has been presented in the context of the space-time approach. The optical 
approach reopens the question of interpretation and reminds us of the effect of pressure on the 
redshift of solar spectral lines. This means that, apart from the gravitational potential, any 
change in mechanical pressure, density and temperature of the medium would give rise to an 
effect of the same kind on the redshift of spectral lines. 

In order to complete the present description, it is necessary to consider the hydrodynamic 
equation. The hydrodynamic equation is 



d_ 

It 



>v) = V( P 4>) - VP + F V . 



(3) 



In addition to gravitational force, we have included pressure term and viscous force. In case of 
a conducting medium with a magnetic field, it is necessary to include the magnetic force term 
in the hydrodynamic equation, leading to the magnetohydrodynamic equation. We restrict 
our discussion to one-dimensional case along the direction of gravitational field. If we neglect 
viscous effects, the hydrodynamic equation takes the form 



d , , d 
(H = 



dt 



dr 



(4) 



Since the time required to propagate a path dr is dr/c to a first approximation, the integration 
of (4) gives 

for the relative change of velocity in the medium which light experiences during propagation 
along the path. The velocity of light at the point of observation thereby becomes 
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as compared with its velocity c at the moment of emission. By Doppler's principle, it can be 
written in terms of frequency as 



f'(r) = f 



1 • ^tI 



AP 
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Optical approach to the redshift effect sheds additional light on its relation to property of 
the medium of propagation. According to optical approach, the redshift effect is a result of 
the velocity of light as affected by pressure gradient including gravitational potential in the 
medium of propagation. From the point of view of electromagnetic waves, light wave in (6) can 
be thought of as a wave in a medium with an index of refraction given by c'(r) = c/n(r). Thus 



n(r) 
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This consideration illustrates how the present picture of redshift offers a natural connection 
with the framework of optics. It has a consequence which is of fundamental importance for 
describing the deflection of light and the radar echo delay from the point of view of optics. To 
reconcile optics and general relativity, however, a vacuum must be understood to exclude even 
a gravitational potential. 

Terrestrial measurements are usually made with respect to a coordinate system fixed in 
the Earth, which rotates uniformly with a constant angular velocity u) relative to the inertial 
system. To an observer in the rotating system, it therefore appears as if the medium is moving 
under the influence of an effective acceleration of gravity^ 4 ) 

g e = g-2(wxv)-wx(wxr). (9) 

The apparent gravitational force acting on the medium is the sum of the actual gravitational 
force, the Coriolis force and the centrifugal force. After this consideration, we must replace <j) 
by an effective potential 4> e with g e in a rotating system. 

The redshift effect was qualitatively in agreement with astronomical observations both in the 
case of the Sun and in the case of white dwarf star like Sirius B where the effect is about thirty 
times larger. However, the quantitative agreement was not very good. While the frequency 
shift in (1) is independent of the point of observation on the solar disk, observations^ 5 ^ have 
shown that the wavelength of spectral lines increases as the point of observation moves toward 
the limb. Furthermore, the solar lines observed at the limb arc definitely asymmetric, having 
pronounced red flanks. There seems to be a systematic change in profile as one approaches the 
limb. In atomic spectra/ 6 ) the broadening of a spectrum line due to pressure has shown that the 
spectrum line observed is spread out more on the long wavelength side than it is on the short. 
With increasing pressure, the mean collision time increases and the time between collisions 
decreases with the result that, as the line is shifted to the red, it is broadened asymmetrically. 
From this point of view, the asymmetry observed in limb lines seems to be of pressure character. 
In fact, Blamont and Roddier^ 7 -* found a complete interpretation of their experimental value at 
the limb when they added to the gravitational redshift the pressure redshift of the Lindholm 
effect. Together with asymmetric profile, their interpretation reminds us of the effect of pressure 
on the redshift of solar spectral lines. 

The controlled experiments using Mdssbauer effect are able to test the gravitational redshift 
to an excellent accuracy. In the experiments/ 8 -* 7— rays in a nuclear resonance passed through 
an evacuated tube or a tube filled with helium along the lines of force of the gravitational 
field. The observed redshifts were shown to be in agreement with the predicted shift in (1) by 
taking account of the effects due to temperature difference between the source and absorber. 
A measurement of redshifts in a rapidly rotating system was shown to fit the dependence of 
g e on w in (9). 

3 Light Bending near the Sun 

The general theory of relativity states that light rays propagated across a gravitational field 
undergo deflection. The theoretical value for the deflection of light rays that just graze the 
Sun's surface is 1.75". The deflection angle is classically measured by comparing the apparent 
positions of stars that happen to lie near the solar disk during an eclipse, when their light comes 
close to the Sun and yet may be detected, with their positions at night six months earlier, when 
these stars lie on opposite sides of the Earth from the Sun, so their light does not pass close 
to the Sun on its way to us. In eclipse of 1919, about a dozen stars in all were observed, and 
yielded values 1.98±0.12" and 1.61±0.31", in substantial agreement with Einstein's prediction. 
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It is perhaps this dramatic result more than any other success that brought general relativity 
to the attention of the general public. 

The Schwarzschild metric, appropriate for the region exterior to a spherically symmetric 
distribution of mass M, is given in the standard form as 

c 2 dT 2 = ^ _ ^y^dt 2 -(l- 1 dr 2 -r 2 (d8 2 + sin 2 9d V 2 ). (10) 

In what follows, we use for the components of the metric tensor the expressions goo( r ) = 
l/g rr (r) = 1 — 2GM/c 2 r. Assuming that the whole motion takes place in the plane ip = 0, 
we obtain as the equations of motion three differential equations. For light rays propagating 
along the geodesic lines, we replace the parameter r by a parameter s describing trajectory. In 
particular, Weinberg^ 10 ) chooses to normalize s so that 

cdt .cdt 
gao~r = 1 111 ffoo^- = constant. (11) 
ds ds 

On the normalization condition the three differential equations can be combined into one dif- 
ferential equation, which is of the same structure as (10). According to Weinberg, the change 
in 9 as r decreases from infinity to its minimum value ro is given by 
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This integral can be evaluated by expanding in the small parameters GM/c 2 r and GM/c 2 r^ to 
first order, giving 4GM/c 2 ro for a light ray deflected by the Sun. 

We now turn our attention to optics, and seek the optical form of expression for the deflection 
of light by the Sun. The gaseous layers surrounding the Sun, through which light propagates, 
are media of spherically symmetric varying refractive index. The deflection of light by the Sun 
can thus be thought as a result of the refraction of light in such gaseous layers from the optical 
point of view. If this thought is reasonable, a similar conclusion will also be reached from this 
point of view. When we check the propagation of rays to this end, we realize that the equation 
for rays in geometrical optics has previously had the form of (12). 

In a homogeneous medium, the refractive index is a constant and the light rays have the 
form of straight lines. Let us consider rays in a medium which has spherical symmetry, i.e. 
where the refractive index depends only on the distance r from a fixed point O: n = n(r). This 
case is approximately realized by the Earth's atmosphere, when the curvature of the Earth is 
taken into account. The light rays are then plane curves, situated in a plane through the origin, 
and along each ray satisfy^ 11 ) 

n(r)r sin ip = constant, (13) 

where ip is the angle between the position vector r and the tangent at the point r on the ray. 
Since V = 7r /2 at the point r of closest approach of the ray to the origin, (13) may also 
be written as n(ro)ro = constant. This relation is sometimes called the formula of Bouguer in 
geometrical optics. If (r, 9) are the polar coordinates of a plane curve, then the angle tp between 
the radius vector to a point r on the curve and the tangent at r is given by 

s '^= ( , r2+rW)1/2 . (14) 

From (13), (14), and Bouguer's formula, the equation of rays in a medium with spherical 
symmetry has been written in the form 
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- 1/2 dr 
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At first sight, we can see a formal equivalence between the equation for rays in geometrical 
optics and the geodesic equation expressed in (12) by Weinberg. The equation for rays is lacking 
a term arising from the difference in path length. For lack of the term, the equation of rays 
corresponds to the case which is obtained when the curvature of the physical space in a region 
of strong gravitational potential is neglected. According to general relativity, the deflection of 
light is due partly to the varying velocity of light and partly to the non-Euclidean character of 
the spatial geometry. Since these are known to contribute equally to the deflection/ 12 ) it can 
therefore be stated that the equation of rays will give a deflection of only half of the correct 
value. 

This result is to be expected on optical grounds, because the non-Euclidean character of 
the spatial geometry has been neglected in optics. In order to compensate for the change in 
light path due to gravitational potential, one can use the notion of optical path. The optical 
path represents the distance light travels in a vacuum in the same time it travels a distance 
in the medium. If a light ray travels in a medium with spherical symmetry, the optical path 
is given by integral over n(r)dr. This means that the radial interval of integration must be 
corrected by multiplication with n(r) to take into account the difference in path length due to 
gravitational potential. Upon integration over n(r)dr instead of the original integration over 
dr, it would yield a result in which the effect arising from the difference in path length is taken 
into consideration. Using the optical path to correct the change in light path, the equation of 
rays is modified to 



The modified equation of rays is in complete agreement in form with (12). From the proposition 
which has just been proved, one may picture what it is to be a curved space in a region of strong 
gravitational potential. When viewed from the present point, the curvature of the physical 
space in the gravitational field of the Sun can best be understood in terms of the medium with 
spherical symmetry in which the path of rays is to be curved. 

A comparison of (16) with (12) identifies n 2 (r) with g rr (r). In the equation of rays, n 2 (r) 
plays exactly the same role g rr (r) has played in the geodesic equation of general relativity. This 
suggests introducing an optical metric tensor n 2 (r) consisting of the gravitational potential plus 
the pressure gradient. Taking only the leading gravitational potential, both of these equations 
give the same result for the deflection of light by the Sun. The optical metric is an extension 
of the Schwarzschild metric to the pressure term, satisfying the eikonal equation as its line 
element. However, the optical metric is based on the interpretation of the deflection of light as 
a result of refraction in the gaseous layers surrounding the Sun. The emphasis should be on 
the fact that there is an agreement in the explicit form between the geodesic equation and the 
equation of rays. 

If we further include a frequency-dependent dielectric constant e{u) defined by Maxwell's 
equation e = n 2 , we can obtain the explicit and integrated form of the optical metric and exhibit 
completely its frequency dependence. Identifying n 2 (r) with the static dielectric constant, wc 
find far above the highest resonant frequency 



where m and e are the mass and charge on the electron, and N is the total number of electrons 
per unit volume. Since the characteristics of the propagation obviously depend on the index 
of refraction n(u>), it is very natural to expect the frequency dependence of the deflection 
so discussed. In fact, Muhleman, Ekcrs, and Fomalont' 13 ) analyzed their experimental data 
by using geometrical-optics techniques in a spherically symmetric refracting medium of index 




(16) 




(17) 
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n(r,w) = 1 + 2GM/c 2 r — 2ne 2 N(r)/muj 2 , where N(r) is the electron-density profile in the 
corona and interplanetary medium. We may draw a parallel between their technique and the 
present approach. 



4 Radar Echo Delay 

In 1964, Shapiro^ 14 -* proposed a fourth test of general relativity. The test involves measuring 
the time delays between transmission of radar signals from Earth to either Mercury or Venus 
and detection of the echoes. Because, according to the general theory of relativity, the speed 
of propagation of light depends on the strength of gravitational potential along its path, the 
time delays are maximum when the inner planets are at superior conjunction and the radar 
signals just graze the solar limb. The maximum round-trip excess time delays are estimated to 
be about 200 ^sec. Such a change, equivalent to 60 km in distance, could be measured over the 
required path length with modern radar equipment by Shapiro and his collaborators/ 15 ^ The 
most reliable of the measured data agree, on the average, with this excess delay predictions of 
general relativity to well within the experimental uncertainty of ±20%. 

According to Weinberg, the time required for light to go from r to r' is given by 
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(18) 



The integral can be evaluated by expanding in the small parameters GM/c 2 r and GM/c 2 r n to 
first order, giving 240 /zsec for the maximum round-trip excess time delay. 

We are now in a position to derive an optical form of expression for the excess time delay. 
Its explicit form will follow from an equation which specifies the path of rays. An accurate 
expression we seek can then be obtained by converting the equation of rays into an equation 
for the path of rays. 

A procedure starts from (14). Substitution of (14) into (13) gives 



n(r)r 2 d9 



(dr 2 +rW)V2 



= constant. 



(19) 



Since the path of rays is ds = (dr 2 +r 2 dd 2 ) 1 / 2 in the polar coordinates of plane curve, this may 
also be written as 

<r)r 2 d9 _ (2Q) 



Solving for ds, we have 



ds 



ds 



constant. 



n(r)r 2 d9 



(21) 



n(r )r 

Bouguer's formula has been used in the above equation. By making use of the integral in (15), 
the variable of integration can be changed from dO to dr, thereby obtaining the result: 
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n 2 (r)r 2 

Hence, by dividing ds by c'(r) — c/n(r), the time of propagation of rays is found to be 
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where c'(r) is the speed of propagation of light in a region of gravitational potential. Although 
the details are altered by the new form of expression, the optical characteristics of (23) remain 
the same as in (15). Thus, for the correct calculation of excess time delay, we must consider in 
addition the difference in the light path of rays due to gravitational potential. 

As discussed in the case of deflection, this requires integrating the resulting equation along 
the optical path. However, it draws a clear distinction between geometrical optics and general 
relativity. This is because (23) has already manifested the form of an integral over optical path. 
The geodesic equation of general relativity teaches us to make correction not only in the velocity 
of light but also in its path length. If wc make a correction in the radial component of the path 
of rays, that is, in the component of the path along the lines of force of the gravitational field, 
the integral in (23) becomes 

2 " 2 ^ r , (24) 
c 

The modified integral is in complete agreement with the geodesic equation formulated by Wein- 
berg. It becomes evident that the equation for rays in geometrical optics also predicts the radar 
echo delay in exactly the same form as given by the geodesic equation of general relativity. There 
is indeed no difference, to first order in GM/c 2 r and GM/c 2 r , between (18) and (24). Again, 
we identify n 2 (r) with g rr (r) in their roles, leading to consider n 2 (r) as an optical metric tensor. 
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5 Plasma Effect of Corona 



In radio astronomy, it is possible to measure the deflection of radio signals by the Sun with 
potentially far greater accuracy than is possible in optical astronomy. At radio frequencies, 
however, it is necessary to analyze the data in terms of a model, in which part of deflection arises 
from general relativity, and the rest is produced by the corona. No prediction can be drawn 
from general relativity since the plasma effect is frequency dependent and the gravitational 
effect is not. In contrast with general relativity, the optical approach affords a direct way to 
calculate the plasma effect of corona on the deflection of radio waves by the Sun. It forms a 
clear view of this optical analogy of general relativity. We are going to evaluate the plasma 
effect of corona. 

As a first important example, we consider the deflection of light as a combination of the 
general relativistic effect and of refraction in the coronal plasma. The expected angular deflec- 
tion can be accurately computed using the frequency-dependent refractive index expressed in 
(17) in the equation of rays (16): 
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(25) 



In order to evaluate this integral, we use in the integrand expansions in the small parameters. It 
is both easier and more instructive to evaluate the integral after the expansions. The integration 
of (25) can be carried to first order in the small parameters with high accuracy. 

The argument of the square root in (25) can be expanded to first order in the small param- 
eters as 
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so (25) gives 



A9 



3- 



-1/2 



1- 



GM 



■ + 



GMr 2ire 2 N{r) 2ire 2 r 2 (N(r) - N(r )) 



c 2 r (r + r ) muj 2 



mui 2 {r 2 — r 2 ) 



(27) 



Consequently, the deflections from the individual effects are combined linearly. Refraction effect 
in the solar corona is now represented by 
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This must be an addition to the general relativity deflection. 

In the Allen-Baumbach model/ 16 ) the electron-density profile in the corona is assumed to 
have the form N(r) — 1.55 x 10 8 (i?o/r) 6 electron/cm 3 , where Rq denotes the radius of the 
Sun. Using more recent results on the corona, Erickson^ 17 ) found that N(r) = 5 x 10 5 (R Q /r) 2 
electron/cm 3 represents the data reasonably well from 4i? Q to 20i? Q . Refraction effect is 
significant where r < 3R & , at which the (i? Q /r) 6 term dominates. Hence, we use the electron 
distribution of the Allen-Baumbach model, resulting in 
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The integration for 59 c is straightforward, and gives 
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(30) 



where cos 9 = r /r. 

The total change in 9 as r decreases from infinity to its minimum value ro and then increases 
again to infinity is just twice its change from oo to ro, that is, 2A9. Hence, the deflection of 
the path of rays from a straight line is given by 59 = 2A9 — ir, which is calculated positively 
if concave toward the Sun and negatively if convex. Putting in the numerical factors, the total 
deflection is 

R. \ 6.24 > 10 1 ' ( 10571-^/?. 



1.75' 
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V 48 
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(31) 



Equation (31) describes interesting behavior of the radiation bending near the Sun. The 
first term represents the general relativistic effect by which the path of rays is bent toward 
the Sun. The second term represents the coronal refraction by which the path of rays is bent 
away from the Sun to the contrary. This is not surprising when we see the difference in sign 
between these terms. Actually, experimental values at radio frequencies of the general relativity 
deflection were determined by fitting, by the method of least squares, the measured data to 
curve of a model bearing difference of sign between these effects. At optical frequencies, coronal 
refraction is extremely small, so it can be neglected. At radio frequencies, however, it plays an 
important part in the deflection, as can be seen when we illustrate the deflection angle as a 
function of frequency for the distances in solar radii of the ray's point of closest approach to the 
Sun's center. Figure 1 shows curves calculated from (31) of the deflection angle as a function 
of frequency for tq/Rq — 1.5, 2, 3, 4. 

The question might be raised as to whether varying velocity of light in the coronal plasma 
also gives rise to a change in path length therein. If we assume that varying velocity of light in 
the coronal plasma does not give rise to a change in path length of rays therein, the radial interval 
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of integration must still be corrected by multiplication with n(r) even in the coronal plasma, 
not with n(r, u) as used in (25). We must then drop the fourth term in the integrand of integral 
in (27), that is, the first term in (28). Coronal refraction thus obtained will be exactly the same 
as what one finds by evaluating the original equation of rays (15) on purely optical grounds. 
Note that there is a complete agreement in the form of expression for the plasma effect between 
(15) and (25) with such an assumption. In fact, the evaluation of coronal refraction from the 
equation of rays (15) was carried out to first order by Bracewell, Eshelman, and Hollweg/ 18 ) 
Their calculation gives 82(i? /r o ) 6 sec for the angular deviation of a ray of frequency 9.6 
GHz in the corona assuming the Allen-Baumbach model. When Erickson's coronal model is 
instead assumed, the angular deviation is given by 0.14(i? Q /r ) 2 sec. Seielstad, Sramek, and 
Weiler^ 19 ) used in data analysis these values as parameters describing refraction effects in the 
solar corona, when they measured the deflection of 9.602 GHz radiation from 3C279 in the solar 
gravitational field using an interferometer at the Owens Valley Radio Observatory. The results 
of their calculation are indeed in exact agreement with what we would obtain from each model 
if we excluded the first term from the integrand of integral in (28) under such a consideration. 
However, on the assumption that varying velocity of light in the coronal plasma also gives rise 
to a change in its path length, we obtain from (31) coronal refraction of 96 (i? Q /r ) 6 sec for 
9.6 GHz frequency. If we used Erickson's coronal model, we would obtain coronal refraction of 
O.21(i? /r o ) 2 sec. The difference of calculation can be written in terms of the metric tensor. 
We have used in (31) the metric tensor of the components 

•9oo = , \, \ and g rr = e(uj)n 2 (r). (32) 

As viewed from the present approach, their calculation corresponds to the case which is obtained 
when the components of the metric tensor are 

■9oo = i \, \ but 9rr = 1 or n 2 (r). (33) 

The reason for this difference is readily understood by referring to the equations of rays (16) 
and (15) from which angular deviations were calculated respectively. 

As a second example of the plasma effect, let us calculate the plasma effect of corona on the 
radar echo delay. The time of propagation of rays is given by (24). To evaluate the dependence 
of the radar echo delay on the frequency being propagated, we must consider (24) with the 
frequency-dependent refractive index in (17): 



t(r ,r') 



n 2 (r ,u!)r 2 
n 2 (r, tv)r 2 



~ 1/2 n 2 (r,uj)dr 



In order to evaluate this integral, we once again use in the integrand the expansions in the 
small parameters to first order. Proceeding in exactly the same way as for (25), (34) gives 



' 1/2 dr r | 2GM | GMr 47re 2 A(r) 2vre 2 r 2 (N(r ) - N (r)) 



c 2 r(r + ro) miv 2 muj 2 (r 2 — r 2 ) 



(35) 

The time required for radar signals to travel to Mercury and be reflected back to Earth is 
2[t(rE,ro) + i(ro, ?"m)]) where te and tm are astronomical radii of the orbits of the Earth 
and the Mercury around the Sun. The round-trip excess time delay is then given by St = 
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2[i(r£7,r ) + t(r ,r M ) - T(r E ,r ) - T(r ,r M )}, where T(r E ,r ) and T(r ,r M ) are the times 
required for radar signals to travel the paths in straight lines at speed c. The distance r of 
closest approach of the radar wave to the center of the Sun is much smaller than the distances 
te and tm of the Earth and Mercury from the Sun. 

Assuming the electron distribution of the Allen-Baumbach model, the integral yields 



AGM 
St ~ — — 



If, instead, we use Erickson's coronal model, we then have 



AGM 
St ~ — = — 



1 + ln^ 



( Ar E r M \ 



2.01 x 10 13 /67rr \ f R e x 2 



P V c / V r o 



(36) 



(37) 



As in the case of deflection, because of the sign difference, the plasma effect of the solar corona 
on the time delay is opposite to what is usually expected from the general relativistic effect. For 
cither equation, the positive terms on the right describe a general relativistic delay in the time 
it takes a radar signal to travel to Mercury and back. In contrast, the negative terms describe a 
coronal plasma contraction, that is, a radar time contraction. Figure 2 shows calculated curves 
of the excess time delay as a function of frequency for ro/i?© = 1.5, 2, 3, 4. 

We now compare Shapiro's calculation with the results obtained from the preceding equa- 
tions. At the beginning of measurement, Shapiro estimated St ~ 140 — 370 ^isec for observations 
of Mercury near superior conjunction with ro ~ ARq at 430 MHz frequency of the Arccibo Iono- 
spheric Observatory. This difference in time delays between the general relativistic effect and 
the coronal plasma effect was nowhere large enough and positive for a really reliable result 
to be obtained solely from Arecibo data. He was thus tried to reduce the plasma effect by a 
factor of almost 400 by using measurements made at 8350 MHz frequency of Haystack radar 
of Lincoln Laboratory. For observation of Mercury with r = ARq at 430 MHz radiation, (36) 
gives St ~ 180 — 126/isec, and (37) gives St ~ 180 — 119/isec. The values obtained for the 
plasma effect are about three times larger than Shapiro's estimate. This is because we have 
made a correction in the radial interval of the path and integrated the resulting equation along 
the optical path bending near the Sun, unlike Shapiro's calculation along the optical path in 
straight line without any corrections. If we assume that the difference in path length is due 
solely to differing gravitational potential, the radial interval of integration must still be cor- 
rected by multiplication with n(r) even in the coronal plasma, not with n(r,u>) as used in (34). 
The excess time delays are then given by 180 — 840 yusec and 180 — 790 fisec, respectively. The 
values so obtained for the plasma effect are just what we should expect if the excess time delays 
were calculated from the equation of rays without any correction on purely optical grounds. 



6 Conclusion 

The optical version of gravitational redshift explains the redshift effect as a result of the velocity 
of light as affected by properties of the medium of propagation. According to this explanation, 
not only gravitational potential but any change in pressure, density and temperature of the 
medium would give rise to an effect of the same kind on spectral lines. Actually, the shifts 
observed in experiments find complete interpretations when the pressure or temperature effects 
are taken into account. From the optical point of view it follows that the effect of pressure 
gradient including gravitational potential can equally be expressed in terms of the medium of 
continuously varying refractive index. To reconcile optics and general relativity, a vacuum must 
be understood to exclude even a gravitational potential. Einstein's interpretation of redshift as 
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a manifestation of time dilatation in a gravitational potential is purely conceptional, which is 
by no means observational. 

It has shown that the equation of rays can also predict the correct values for the deflection of 
light and the radar echo delay by the Sun. To be reconciled with general relativity, I have made 
a correction to the path length in the equation for rays in geometrical optics. An agreement 
in form between the equations of rays and the geodesic equations suggests introducing an 
optical metric of n 2 (r) consisting of the gravitational potential and the pressure gradient in 
the medium of propagation. These results are obtained in the explicit form on a comparison 
of the equations of rays with the geodesic equations formulated by Weinberg. In the light of 
this fact, Weinberg's formalism has opened a door to introduce a relationship between general 
relativity and geometrical optics. Indeed, it was the formal agreement between (12) and (15) 
that enabled the present approach to be proposed. 

The equation of rays provides a theoretical curve for observation that any beam of radiation 
is deflected during its passage near the Sun as a result of the gravitational effect and of refraction 
in the coronal electron plasma. The optical analogy of general relativity may be supposed to 
serve as a means of describing additional gravitational effects involving optical phenomena. But 
it should be noted that this analogy is made to explain the general relativistic effects from the 
point of view of optics. In a comparison with the previous estimates, special attention is given 
to a change in the path length of rays. As in the case of gravitational effect, I have assumed in 
the path length a change caused by the plasma effect. 
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Note added in proof 



I should like to express gratitude to a reviewer for corrections, criticism, and 
language corrections of the early draft of this paper. I guess the reviewer to be 
Prof. M. D. Kruskal. The reviewer's comments are so instructive that I would like 
to accompany his reports with this paper. 

[1] The rcdshift equation, Eq. 1 should hold for any acceleration. If derivable 
from a potential one could postulate that gh/c be replaced by the right member of 
Eq. 5 and then define n(r) by Eq. 8. This would avoid the potential confusion of 
fluid drag effects, as in Fizeau's experiments, and also divorce the proposed effect 
from fluid flow per se. As the effects all stem from relative accelerations of reference 
frames, a clear description of what constitutes proper frames within and outside a 
"fluid" is needed. 

[2] I do not think that it is always possible to attribute the redshift to a relative 
velocity between an inertial frame and a frame moving with a fluid. It may be 
possible that a relative acceleration between frames causes a local redshift and a 
cumulative time delay. I think that an appeal to the principle of equivalence is still 
necessary in order to include all causes of acceleration as causes of redshift. 

[3] Without doing all of the manipulations required for proof, I suspect that 
one could retrieve the redshift implied by Eq. 5 from the field equations of general 
relativity. The terms beyond gravitational potential should be contained within the 
stress-energy tensor of the right member of the field equations. This is no criticism 
of Eq. 5; in fact it is good to have an intuitive way of arriving at the redshift result. 

[4] The author describes Einstein's interpretation of redshift as a manifestation 
of time dilation in a gravitational potential as "rather unusual" . I believe that this 
is unwarranted. At this point in the development of the author's optical analogy he 
has tied everything to changes caused by properties of the "medium" of propagation. 
Since there are effects clearly due to a gravitational field in the absence of any other 
"medium", and other interpretations are lacking at this point, I believe that it is 
incorrect to call Einstein's interpretation "rather unusual" . 

Reply: I learn much from his comments. But I cannot agree to the point [4]. 
It reveals a difference of standpoint looking at the relativity theory between the 
current paradigm (reviewer) and the opponent (author). What the experiment has 
actually demonstrated is a change in wavelength or frequency with gravitational 
potential, not a change in rate of clock with it. One can see from the original paper 
that Einstein presented the interpretation of the rcdshift effect as a result of fitting 
the predicted speed of light d = c(l + gh/c 2 ) even deliberately into the postulate 
of the constancy of the speed of light. I think that Einstein put the cart before the 
horse in its interpretation, apart from the controversy of whether the postulate is 
consistent or inconsistent. 
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